Chaotic instantons and enhancement of tunneling in kicked double-well system with 

time-reversal symmetry 
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Chaotic instanton approach is used to describe dynamical tunneling in kicked double well sys- 
tem. Effective Hamiltonian for the kicked system is obtained using matrix expansion formula for 
operator exponent and exploited to construct an approximation for chaotic instanton solution. This 
approximation is used for derivation of the ground quasienergy splitting dependence on both the 
perturbation strength and frequency. Results of numerical calculations for corresponding ground 
quasienergy splitting dependencies based on Floquet theory are in good agreement with the derived 
analytical formula in a wide range of perturbation parameters. 
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I. INTRODUCTION 

Investigation of the influence of small perturbation on 
the behavior of the nonlinear dynamical systems attracts 
permanent interest for the several last decades [H-ll] . The 
connection between the semiclassical properties of per- 
turbed nonlinear systems and purely quantum processes 
such as tunneling is a reach rapidly developing field of 
research nowadays [3, 5]. Our insight in some novel phe- 
nomena in this field was extended during the last decades. 
The most intriguing among them are the chaos assisted 
tunneling (CAT) and the closely related coherent de- 
struction of tunneling (CDT). 

The former in particular is an enhancement of tunnel- 
ing in the perturbed low-dimensional systems at small 
external field strengths and driving frequencies [B-Q- 
This phenomenon takes place when levels of the regu- 
lar doublet undergo an avoided crossing with the chaotic 
state [Hi H3 • At the semiclassical level of description one 
considers tunneling between KAM-tori embedded into 
the "chaotic sea". The region of chaotic motion affects 
tunneling rate because compared to direct tunneling be- 
tween tori it is easier for the system to penetrate pri- 
marily into the chaotic region, to travel then along some 
classically allowed path and finally to tunnel onto an- 
other KAM-torus [SHI- The later, CDT phenomenon, 
is a suppression of tunneling which occurs due to the ex- 
act crossing of two states with different symmetries from 
the tunneling doublet [14]] • In this case tunneling time 
diverges which means the total localization of quantum 
state on the initial torus. 

CAT phenomenon as well as CDT were experimen- 
tally observed in a number of real physical systems. The 
CAT observation between whispering gallery-type modes 
of microwave cavity having the form of the annular bil- 
liard was reported in Ref. [15j. The same phenomenon 
for ultracold atoms was experimentally investigated in 
Refs. [Ill I!?]]- The study of the dielectric microcavities 
provided evidences for CAT as well [18|. Both CAT and 
CDT phenomena were observed in two coupled optical 
waveguides [HI, [2(| ■ Recently experimental evidence of 



coherent control of single particle tunneling in strongly 
driven double well potential was reported in Ref. (2l| . 

A number of effective approaches [22rf24| are used for 
analysis of the classical and quantum properties of the 
perturbed systems [25|, [26[. The most common meth- 
ods which are used to investigate the interplay between 
semiclassical properties of perturbed nonlinear systems 
and quantum processes are numerical methods based on 
Floquet theory [|, [H, H3 and Random Matrix Theory 
[28j . Among other approaches we would like to men- 
tion the scattering approach for billiard systems [29|, H3| 
and approach based upon the presence of a conspicuous 
tree structure hidden in a complicated set of tunneling 
branches [3TM33I ] - 

In this paper we will consider the original analytical 
approach based on instanton technique. Enhancement 
of tunneling in system with external force in framework 
of this approach occurs due chaotic instantons which ap- 
pear in perturbed case. This approach was proposed in 
p3 - l36| and used in [37j ■ Chaotic instanton approach 
will be developed further here using effective model and 
exploited for description of the enhancement of tunnel- 
ing in the kicked double well system. The main purpose 
of the present study is to prove the ability of developed 
chaotic instanton approach to give quantitative analyti- 
cal description of tunneling well agreed with independent 
numerical calculations based on Floquet theory. It will 
give additional support and pulse for the further devel- 
opment of analytical methods to investigate tunneling 
phenomenon in quantum systems with mixed classical 
dynamics. Alternative approach based on quantum in- 
stantons which are defined using an introduced notion of 
quantum action was suggested in (38l . [39| . Analytical ap- 
proach to describe tunneling in perturbed systems based 
on nonlinear resonances consideration was developed in 
[25l Eoj . A theory for dynamical tunneling process us- 
ing fictitious integrable system was recently given in Ref. 

El- 
Double well potential is a special toy model among 
other nonlinear systems. This system is a simple and 
nonlinear at the same time. It is convenient to use this 
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model for tunneling analysis. This system is well studied 
in the nonperturbed case, e.g. on the base of instanton 
technique [H, El] or WKB method [HJ - Double well po- 
tential is often used for description of processes which 
occurred in wide range of real ph ysic al systems: such as 
flipping of the ammonia molecule [45j , transfer of protons 
along hydrogen bonds in benzoic-acid dimers at low tem- 
peratures [46l . |47| and macroscopic quantum coherence 
phenomena in superconducting quantum interference de- 
vices [48l - l50j and nanomagnets (SJ HH . Perturbation in 
this paper is regarded in the form of the periodic kicks. 
One of the attractive features of this type of perturba- 
tion is the extensively-investigated simple quantum map 
which stroboscobically evolves the system from kick n to 
kick n + 1. Kicked systems are recently used for experi- 
mental realization of a such novel concept as a quantum 
ratchet pH^. 

The main aim of this study is a development of the 
chaotic instanton approach on the base of the construc- 
tion of model system effectively approximating original 
chaotic dynamics by some averaged regular trajectories. 
Tunneling in kicked double well system is investigated us- 
ing developed chaotic instanton approach. For this pur- 
pose nonperturbed trajectories of the system in imagi- 
nary (Euclidean) time are analyzed in section HU Ef- 
fective model approximating euclidean chaotic dynamics 
of the original system is constructed using matrix ex- 
pansion formula for one period evolution operator in sec- 
tion [IlTJ Chaotic instanton and its approximation in con- 
structed time-independent effective model are discussed 
in section IIV1 Results obtained by means of the effec- 
tive model are used in section [V] to derive analytical for- 
mula for lowest quasienergy doublet splitting dependence 
on perturbation parameters. Numerical calculations are 
performed to check the validity of this formula in the 
section IVU 



II. KICKED DOUBLE- WELL POTENTIAL AND 
NONPERTURBED TRAJECTORY ACTION 

Hamiltonian of the particle in the double- well potential 
can be written down in the following form: 



Hn = 



P 

2m 



4 2 

ao x — a,2 x , 



(1) 



where m - mass of the particle, ao, a2 - parameters of the 
potential. We consider the perturbation of the kick-type 
and choose it as follows: 



£ X 



2 E 



S(t-nT), 



(2) 



where e and T are perturbation strength and period, re- 
spectively, t - time. Dependence of the perturbation on 
coordinate was chosen in the form of x 2 in order to pre- 
serve spatial symmetry in the perturbed system. Hamil- 
tonian of the perturbed system is the following: 



H = Hn 



(3) 



We use Euclidean time in this section to work with 
instantons in section IIVI Now we implement Wick rota- 
tion (t — > —it) and define Euclidean Hamiltonian in the 
following form: 



ft — rtn 



e x 



E 

n= — oo 



S(t-uT), 



(4) 



where T~Lq - nonperturbed Euclidean Hamiltonian which 
can be written down as follows: 



njE P 4,2 

Hn = a x + a 2 x . 

2m 



(5) 



It is seen that all perturbed trajectories, as well as 
chaotic instantons, are composed of nonperturbed tra- 
jectories' parts connected by jumps which are caused by 
kicks. Thus primarily we will investigate nonperturbed 
dynamics in this section and obtain the expression for 
nonperturbed trajectory action near the separatrix. We 
use the action-angle variables defined in the case of one- 
dimensional systems. Using standard technique [l], sec- 
tion 1.2] for the nonperturbed Hamiltonian ([5]) we obtain 
expressions for dependence of these variables on energy E 
and coordinate x in double well system in Euclidean time 
and write down them in the following form: 



J(E) 



= 2_a^Vm I J^ao 



37TO0 



E 



e(E, x ) 



x \ L{ X )- . l-44£ K( X ) 



-F 



(6) 




1-4^-E 



J 



where J and 9 - action and angle variables, respectively, 
K(x) and L(x) - full elliptic integrals of the first and sec- 
ond kinds, respectively, F (</>, x) - elliptic integral of the 
first kind, where 4> is integral argument, x ~ its modulus. 
Modulus x can be expressed in the following form: 



l-4<a-E 
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We expand expression ([6]) in powers of the (Euclidean) 
energy difference from the separatrix £ (£ = Ei ns t — E, 
Einst = ai/^ao) is the energy on the separatrix) to ob- 
tain analytical ex pre ssion for nonperturbed trajectory ac- 
tion (see [34l. [35l l55j ) . For this purpose we use the rep- 
resentation of full elliptic integrals with close to unity 
modulus x [HI chapter 21.6]: 



L(X) = 1 + \ (A - \ I \ 
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^(x) = A + i(A-l )x ' + A(A-0 x ' 2 
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where A = ln(4/y / x 7 ) and \' — 1 — X- Substituting these 
expressions in (|6]) and neglecting terms higher than linear 
one we obtain the following estimative relation between 
the nonperturbed trajectory action and the Euclidean 
energy difference from the separatrix: 



S[x(t,£)] = irJ(E inst -£) = S[x ins t(r, 0)] - a, — £, 

V a 2 

(7) 

where S[x inst {r, 0)] = <Si„s« = 2^/ma^ 2 /(3 a Q ) - nonper- 
turbed instanton action, a = (1 + 181n2)/6 - numerical 
coefficient. 



III. EFFECTIVE TIME-INDEPENDENT 
MODEL FOR THE KICKED SYSTEM 

Now lets construct the effective time-independent 
Hamiltonian for the double well system with the per- 
turbation of the kick-type. Initially we will calculate this 
Hamiltonian in the real time. Effective Hamiltonian in 
Euclidean time will be obtained at the end of the sec- 
tion and it will be used for the Euclidean phase space 
analysis in the next one. We will construct the effec- 
tive time-independent Hamiltonian for the system under 
investigation using the following definition [22j : 



ff J 



e -iHoT/2 e - 



iex 1 -iH T/2 

& 5 



(8) 



Nonperturbed Hamiltonian (JTJ is the first term. The sec- 
ond term has the same order (zero) as the nonperturbed 
part on small parameters (e, T). Terms having the first 
order are absent in the expression (JSJ . Next terms which 
are shown in the formula have the second order on pertur- 
bation parameters. Terms with higher orders are omitted 
in the last expression. We restrict our consideration by 
two first terms and neglect the others. Thus effective 
Hamiltonian for the kicked double well system describ- 
ing the tunneling dynamics between wells can be written 
down in the following form: 



H, 



eff 



P 

2m 



CIqX 
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C V 

2^ 



(10) 



This Hamiltonian is coincided with nonperturbed Hamil- 
tonian when parameter a 2 is replaced by a 2 . The later 
is defined as follows: 



a 2 (e,v) = a 2 - —. 

Zn 



(11) 



In contrast to the kicked system ^ the effective system 
with Hamiltonian (|10| is time-independent. We should 
introduce one more restriction for the perturbation pa- 
rameters variation. This restriction follows from our as- 
sumption that ordinary instanton approach should be 
valid for the effective potential. Condition for ordinary 
instanton approach applicability is the following |43| : 



where RHS is a one-period evolution operator of the 
kicked system. During the kicks the influence of the 
nonperturbed Hamiltonian is negligible in comparison 
with perturbation. The moment of the kick is chosen 
in the middle of the period in order to conserve the time- 
reversal symmetry in the system (i — > —t) Q. We re- 
strict our consideration by sufficiently small values of 
both the perturbation strength and period (e < 0.1 and 
T < 2-7T /wo, where ujq - oscillation frequency near the 
bottom of the wells). Deriving H e ff from definition (jHJ 
and using matrix expansion formula for the exponents 
we can write down effective Hamiltonian in the following 
way: 



H, 



eff 



±ln 
T 



1 - iH Q T 



H 2 T 2 



cT 



(H x 2 + x 2 H Q ) + 0(3) 



where we show terms up to the second order on e and T. 
Finally we obtain the effective time-independent Hamil- 
tonian for the kicked system expanding the logarithm in 
the series: 
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Hn 
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We obtain the following restriction in the case considered: 



'-- v ^ (e v) max = 2tt 



0,2 



(12) 



Thus we should use intermediate values of the pertur- 
bation frequency. Frequency should be not small to use 
matrix expansion formula and not larger than (e v) m ax/^- 

Effective potentials for the set of the perturbation 
strength and frequency values are shown in the figure [1] 
The model parameters for all figures in the paper are the 
following m = 1, ao = 1/128, a 2 = 1/4. The two 
lowest energy levels in the effective system are shown 
by horizontal lines in the figure. One can see that for 
parameters' values satisfying the condition (TT21) the two 
lowest levels are quasidegenerate (e = and e = 0.08). 
Wells became closer and barrier height is decreased when 
we increase perturbation parameters. It should enhance 
tunneling rate in kicked system. These consequences will 
be checked in following sections for quasienergy spectrum 
and tunneling process in the system. Two lowest lev- 
els are shifted one from another and potential has one 
well for large values of the perturbation strength and fre- 
quency (e = 0.5 and v = 4 in case shown in the figure [1}. 
These perturbation parameters are not satisfy the condi- 
tion CO 
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FIG. 1: Potential and two lowest levels for the effective model. 
Perturbation frequency v — 4, perturbation strength e = 
(thick solid line), e = 0.08 (thin solid line) and e = 0.5 (dashed 
line). The model parameters are m = 1, do = 1/128, 
a 2 = 1/4. 



Now we will perform the transformation for effective 
Hamiltonian (|10[) into Euclidean time and obtain the fol- 
lowing expression: 



ttE 
U eff 



P 

2m 



a 2 



c v 
2^ 



(13) 



This Hamiltonian will be used in the following section 
to analyze the perturbed system phase space in Euclidean 
time and to construct an approximation for chaotic in- 
stanton solution. 



IV. CHAOTIC INSTANTONS AND EFFECTIVE 
MODEL EUCLIDEAN PHASE SPACE 

Let us regard a possibility to describe properties of the 
classical motion in the kicked double well system in Eu- 
clidean time using effective model (fT3"|) .We will conduct 
numerical simulations of the particle dynamics both in 
the kicked system and in the framework of effective model 
for the same set of initial conditions. The deviation be- 
tween the trajectories in phase space when certain time 
expires will be a criterion of the applicability of the time- 
independent approximation in the classical case. 

We will consider kicked system Q with the perturba- 
tion parameters e = 0.02 and v = 7 and the effective 
model with parameter ai which is defined by the expres- 
sion (jlll) . Separatrix in the effective model is shown in 
the figure [2] by thick solid line while separatrix in nonper- 
turbed system by dashed line. Effective trajectories with 
energies which are less then on the separatrix are shown 
by thin solid lines. Comparison of the classical motion 
on one period of the perturbation in the effective (thick 
solid lines) and kicked (thin solid lines) systems for the 
set of initial conditions (thick points) are shown in the 
inset (a) in the figure Trajectory in the effective model 



is a smooth line in the phase space. The trajectory for 
the perturbed system is a broken line due to the kick in 
the middle of the period. Trajectories for two systems lie 
close to each other as it is shown in the inset (a) in the 
figure [2j Inset (b) shows chaotic trajectory of the kicked 
particle near the turning point of the effective separatrix. 
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FIG. 2: Effective model phase space for the kicked system 
with perturbation parameters e — 0.02 and v — 7. Separatrix 
in the effective model (thick solid line) and in the nonper- 
turbed system (dashed line) are shown in the figure. Com- 
parison of the particle classical motion on one period of the 
perturbation in effective (thick solid lines) and kicked (thin 
solid lines) systems from the set of initial conditions (thick 
points) are shown in the inset (a). Inset (b) shows chaotic 
trajectory of the kicked particle near the turning point of the 
effective separatrix. 

In order to check applicability of the effective model for 
the description of the perturbed system phase space we 
perform numerical simulations of the particle dynamics 
using direct simulations for the kicked system and calcu- 
lations using effective Hamiltonian (fl3)) from the set of 
initial points. Energy range for this set was chosen from 
the bottom of the wells (Eb ott om = 0) to the effective 
model separatrix energy which can be written down in 
the following way: 



E top = a 2 (e, t/) 2 /(4a ). 



(14) 



Fifty energy values were chosen in this energy interval. 
One hundred of initial points were randomly evaluated 
for each energy value. Effective model sufficiently strictly 
describe the classical dynamics of particles in the kicked 
double well system for all particle energy range. Mean 
arithmetical deviation in energy is less then 0.1% of the 
barrier height. Thus phase portrait of the effective model 
coincides with the stroboscopic map for the kicked sys- 
tem besides separatrix region. Narrow stochastic layer 
should be placed near the separatrix in the kicked sys- 
tem. Inset (b) shows this layer near the turning point of 
the effective separatrix. 

Euclidean equations of motion of the particle in the 
nonperturbed double-well potential (e = 0) have a well 
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known solution - instanton. This solution is used for 
calculation of the ground energy splitting in the system 
without perturbation [H, |43| and explains the rate of 
the tunneling process in it. Another solutions of the 
Euclidean equations of motion besides ordinary instan- 
ton are required to explain dynamical tunneling in per- 
turbed system. Perturbation destroys the separatrix and 
some trajectories in its vicinity go to infinity. Narrow 
stochastic layer is formed nearby the nonperturbed sep- 
aratrix due to the perturbation. "Chaotic instanton" is 
appeared in this layer. It consists of nonperturbed tra- 
jectory's parts connected by jumps. Chaotic instanton is 
the closest to the destroyed separatrix perturbed trapped 
trajectory. Thus it plays a dominant role in tunneling in 
perturbed system. This configuration will be investigated 
analytically in this section. 

Chaotic instanton from dynamical point of view is a 
set of nonperturbed trajectory parts connected by jumps. 
These jumps are induced by kicks with period T. Energy 
of the nonperturbed trajectories will be changed after a 
kick. Effective model can be used for analysis of the 
kicked system phase space. Thus we can use effective 
model separatrix as an approximation for the chaotic in- 
stanton in the kicked system. Using this approximation 
we can write down an expression for description of chaotic 
instanton in phase space: 



p s (x) = 4 2m 



Q2(e, v) 2 
4a n 



d 2 x 2 



(15) 



This approximation, i.e. effective separatrix, is shown by 
thick solid line in the figure [2] Energy of the particle 
which is moving over the chaotic instanton trajectory is 
changing from maximum energy when particle momen- 
tum is zero to minimum at point where particle coordi- 
nate is zero. 

Lets write down expressions for these energy borders in 
explicit form. Chaotic instanton energy is maximum for 
Ps(%±) — 0- This point is a turning point for the effective 
model (x± — ± y / & 2 /(2ao) ) ■ Using condition for turning 
points we can write down expression for maximum energy 
in the following form: 



Emax — (P s 
— Einst 



= 0,x±) 
16ao7r 2 



4a 



16ao7r 2 



(16) 



where Ei nst = a|/ (4ao) - energy of nonperturbed instan- 
ton. We restrict our consideration by small values of the 
perturbation strength e and intermediate values of the 
perturbation frequency v. Thus we neglect the second 
term in the expression (|16[) . 

Chaotic instanton minimal energy is calculated from 
the condition x = for expression (I15[) . This energy 
coincides with the energy on the effective model separa- 
trix (1141) . Minimal energy can be written down in the 



following way: 



Emin = HZ( Ps {0),x = 0) = E top = 2 2 (c > i/)7(4ao) 



of _ (e^) 2 
4ao 4ao7r 16ao7r 2 



Ei, 



4an7r ' 



(17) 



where we have neglected terms higher than the linear 
one. 

The nonperturbed trajectory with energy E max is close 
to the nonperturbed separatrix. The nonperturbed tra- 
jectory with energy E m i n coincides with the effective sep- 
aratrix for point x = and p s (0). Using the last two ex- 
pressions (|16p and (fl~7| we obtain the formula for energy 
range of the chaotic instanton: 



inst. E max 



a 2 tv 
4ao 7r 



(18) 



Expression (|18p will be used in the following section 
in order to obtain analytical formula for the lowest 
quasienergy doublet splitting dependence on the pertur- 
bation parameters in the kicked system. 



V. GROUND DOUBLET QUASIENERGY 
SPLITTING FORMULA 

The lowest doublet energy splitting in two loop ap- 
proximation in the nonperturbed double well potential is 
the following ([57j and review [43j]): 



AE Q = 2uq\I ^- \J S inst exp 



'Si; 



71 1 

T2 Sint 



(19) 



where ujq - oscillation frequency near the bottom of the 
wells , Si ns t - nonperturbed instanton action. 

Ground doublet quasienergy splitting (A77) in the 
kicked system in the framework of our approach is 
expressed in terms of chaotic instanton action (S c h) 
through the formula which is similarly to (|19l) : 



A77 




(20) 



where chaotic instanton action can be calculated by aver- 
aging the nonperturbed trajectory action ([7]) over energy 
from minimum to maximum for the chaotic instanton en- 
ergy (E mm and E max , respectively): 



S c h = 



1 



M4 E 
L - yrL ch.inst. 



S(E)dE. 



This integral can be transformed to the integral over the 
energy difference £ = E inst — E and calculated directly. 
As a result we have the following expression: 



S, 



ch 




(21) 
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Now we can write down analytical formula for 
the ground quasienergy levels splitting using expres- 
sions (US), CE]), ([201) and ([21]): 



Ar](e, v) = AE e 



k e v 



(22) 



where 



k 



■ \fm a 2 
8nao 



Tunneling period in the kicked double well potential is 
expresses in terms of ground quasienergy levels splitting 
through the following formula: 



2tt 
At]' 



(23) 



Increasing of the perturbation parameters gives exponen- 
tial rise to ground quasienergy splitting and to the tun- 
neling frequency {v tU n{^^v) — Ar/(e,h 1 )). The last ex- 
ponential factor in the expression (|22|) is responsible for 
the tunneling enhancement in the perturbed system. In 
nonperturbed case formula (|22p coincides with the ex- 
pression (|19l) . Formulas (|2"2")l and (|23l) will be checked in 
numerical calculations in the next section. 



VI. NUMERICAL CALCULATIONS 

For the computational purposes it is convenient to 
choose the eigenvectors of harmonic oscillator as the basis 
vectors. In this representation matrices of the Hamilto- 
nian ([1} and the perturbation ^ are real and symmetric. 
They have the following forms (n > m): 



11 m n — u m n 



hhj ( n H — I + 



9 (Z 



2 V2 



ga (2m 2 + 2m + 1) 



a' 2 (2m + l))] 



S m +2 n ^ (ff a o(2m + 3) - a! 2 ) \J (m + l)(m + 2) 



ao5 



+ S m +4, n - j- Vim + l)(m + 2)(m + 3)(m + 4), 
V m V(m + l)(m + 2) + S m n (2m + 1) ) , 



where g = fi/raw and a 2 = a 2 + mw 2 /2, ft is Planck 
constant which we put equal to 1, uj - frequency of the 
basis harmonic oscillator which is arbitrary, and so may 
be adjusted to optimize the computation. We use the 
value uj = 0.2 with parameters m — 1, <2rj = 1/128, 
a 2 = 1/4 which are chosen in such a way that nonper- 
turbed instanton action is large enough for energy split- 
ting formula for nonperturbed system to be valid and not 
too big in order to decrease errors of numerical calcula- 
tions. The matrix size is chosen to be equal to 200 x 200. 
Calculations with larger matrices give the same results. 
System of computer algebra Mathematica was used for 
numerical calculations. 
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FIG. 3: Quasienergy spectrum for the ten lowest average en- 
ergy levels. Thick lines - doublet with the minimal average 
energy. Perturbation frequency v = 5. 



We calculate eigenvalues of the one-period evolu- 
tion operator (RHS of the expression (|8])) and obtain 
quasienergy levels (jik) which are related with the evo- 
lution operator eigenvalues (Afc) through the expression 
r/k = i InAfe/T. Then we get ten levels with the low- 
est one-period average energy which are calculated using 
the formula (vi\H + V/T\vi) (\vi) are the eigenvectors 
of the one-period evolution operator). The dependence 
of ten lowest levels quasienergies on the strength of the 
perturbation is shown in the figure |31 Quasienergies of 
two levels with the minimal average energy are shown by 
thick lines in the figure |3] They are too close to each 
other to be resolved in the figure due to very small split- 
ting. 

Performed numerical calculations give the dependence 
of the ground quasienergy splitting both on the strength 
(figffla)) and the frequency (figHJb)) of the perturba- 
tion. Results of numerical calculations are plotted in 
the figure 2] by points. Axis Arj is shown in logarithmic 
scale. Obtained dependencies are exponential as it was 
predicted by chaotic instanton approach and obtained 
analytical formula (|22|) . 

Analytical results are plotted in the fig- 
ures |3l (a) and (b) by straight solid lines. Numerical 
points lie close to these lines. The agreement between 
numerical calculations and analytical expression is good 
in the parametric region considered. 

Now lets perform numerical simulations for the tunnel- 
ing process in the kicked double well system and check 
an applicability of the formulas (|2"2")l and (|2"3")l for this 
process. For this purpose we regard the double well po- 
tential ([1]) with the same parameters and basis vectors as 
for previous calculations. We take a symmetric superpo- 
sition of two lowest nonperturbed states as a initial wave 
packet. These packet is localized in left well of poten- 
tial. Numerical simulations we provide by multiplying 
initial wave function by one period evolution operator 
e -iH T/2 e -iex e -iH T/2^ resm j- g of numerical simu- 
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FIG. 4: Quasienergy splitting as a function of the strength (a) 
and frequency (b) of the perturbation. Lines - analytical for- 
mula (I22|l . points - numerical results. 




FIG. 5: Quantum mechanical tunneling in kicked double well 
potential. Perturbation parameters: v = 4, (a)e = 0, (6)e = 
3-1CT 2 . 



lations for the two sets of the perturbation parameters 
are shown in the figure [5] The dependence of the local- 
ization probability of the wave packet on the coordinate 
and time is presented in figures. Minima of the nonper- 
turbed double well potential ((TJ) are situated in points 
x = —4 and x — 4. Tunneling between these points in 
nonperturbed system is demonstrated in the figure [SJa) . 
Evolution of the initial wave packet in perturbed case is 
shown in the figure EJb). Perturbation parameters for 
these simulations are v = 4 and e = 3 • 1CP 2 . They are 
chosen in such a way to speed up a tunneling in two times 
in comparison with nonperturbed system. Figures [5ja) 
and Ob) demonstrate this enhancement. Fourier anal- 
ysis of the dependence of the localization probability of 
the wave packet in left well on time in perturbed case 
confirms analytical assumptions mentioned above. 

In order to check applicability of the developed ap- 
proach we carry out a series of the numerical calculations 
for wide range of the perturbation parameters. Result of 
the analysis is performed in figure [6] Region of quantita- 
tive agreement between analytical and numerical results 
is shown by black color in the figure. There are two re- 
strictions of the developed approach applicability. The 



Chaotic instanton approach applicability error in percents 




0.02 0.04 0.06 0.08 0.1 

Perturbation strength, ^ 



FIG. 6: Comparison error of the analytical formula ()22[) with 
results of numerical calculations in percents. 



first one is that the perturbation frequencies have to be 
not too small. Otherwise the chaotic instanton approach 
is broken due to restriction for value of perturbation pe- 
riod (e < 0.1 and T < 2t:/ujq) which was introduced in 
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section Hill Another restriction for analytical predictions 
is a condition which imply the limit for product of per- 
turbation strength and frequency (see [P2"|) . These two 
restrictions explain accurately the figure [5] 

Inverse sign in the expression of the perturbation ([2]) 
will induce exponential suppression of tunneling in the 
system. It can be demonstrated numerically as well. 

VII. CONCLUSIONS 

Chaotic instanton approach allows to describe analyti- 
cally the influence of the perturbation on quantum prop- 
erties of nonlinear systems. The kicked double well sys- 
tem is regarded as a toy model to compare quantitative 
analytical predictions with the results of numerical cal- 
culations. 

Chaotic instanton is the solution of the Euclidean 
equations of motion of the perturbed system. This con- 
figuration is responsible for the enhancement of tun- 
neling far away from the exact (avoided) level cross- 
ings. Time-independent effective model is used for reg- 
ular approximation of the chaotic instanton solution in 
order to take into account it's contribution to the ground 
quasienergy doublet splitting. The effective Hamiltonian 
is constructed for the kicked double well system using 
matrix expansion formula for operator exponent for the 
small values of the strength and period of the pertur- 



bation. Dynamical simulations for the perturbed sys- 
tem and effective model show that effective Hamiltonian 
can be used for the description of the perturbed system 
Euclidean phase space. The chaotic instanton approxi- 
mation was constructed and exploited to obtain the en- 
ergy range for the chaotic instanton trajectory. Formula 
for ground quasienergy levels splitting was evaluated av- 
eraging trajectory action in the obtained energy range 
in the framework of chaotic instanton approach. This 
formula predicts exponential dependence of the ground 
doublet splitting on both the perturbation strength and 
frequency. 

Numerical calculations for quasienergy levels depen- 
dence on the perturbation parameters and simulations for 
tunneling dynamics are performed to check the validity 
of the obtained analytical formula. Results of numerical 
calculations for the quasienergy spectrum confirm the ex- 
ponential dependence of the ground splitting on both the 
perturbation strength and frequency. They are in good 
agreement with the derived analytical formula (|22[) . Sim- 
ulations of the tunneling dynamics of the wave packet in 
the kicked double well system demonstrate exponential 
tunneling enhancement as well. Applicability of chaotic 
instanton approach was tested in a series of numerical 
calculations. Sufficiently wide range of perturbation pa- 
rameters was found suitable for developed approach ap- 
plication. 
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